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We discuss production of charm quarks, mesons as well as nonphotonic 
electrons in pp scattering at RHIC. The distributions in rapidity and trans- 
verse momentum of charm and bottom quarks/ antiquarks are calculated 
in the fcj-factorization approach. The hadronization of heavy quarks is 
done by means of fenomenological fragmentation functions and semilep- 
tonic decay functions are found by fitting semileptonic decay data. Good 
description of the inclusive data at large transverse momenta of electrons is 
obtained and a missing strength at small transverse momenta of electrons 
is found. 

In addition we discuss kinematical correlations between charged leptons 
from different mechanisms. Reactions initiated by purely QED 7*7*-fusion 
in elastic and inelastic pp collisions as well as diffractive mechanism of 
exclusive cc production are included. A good description of the dilepton 
invariant mass spectrum of the PHENIX collaboration is achieved. Distri- 
butions in the dilepton pair transverse momentum and in azimuthal angle 
between electron and positron are presented. 

A new mechanism of exclusive production of cc is discussed. Corre- 
sponding results are shown and the possibility of its identification is dis- 
cussed. 

We discuss also production of two pairs of cc within a simple formalism 
of double-parton scattering (DPS). Very large cross sections, comparable to 
single-cc production, are predicted for LHC energies. Both total inclusive 
cross section as a function of energy and differential distributions are shown. 
We discuss a perspective how to identify the double scattering contribution. 
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1. Introduction 

In recent years the PHENIX and STAR collaborations have measured 
transverse momentum distribution of so-called nonphotonic electrons [TJ [2] . 
The dominant contribution to the nonphotonic electrons/positrons comes 
from the semileptonic decays of charm and/or beauty mesons. Formally 
such processes can be divided into three subsequent stages. First cc or bb 
quarks are produced. The dominant mechanisms being gluon-gluon fusion 
at higher energies or quark-antiquark annihilation close to the threshold. 
Next the heavy quarks/antiquarks fragment into heavy charmed mesons 
D,D* or B,B*. The vector D* and B* mesons decay strongly producing 
pseudoscalar D and B mesons. Finally the heavy pseudoscalar mesons decay 
semileptonically producing electrons /positrons. 

The hadronization of heavy quarks is usually done with the help of 
phenomenological fragmentation functions with parameters adjusted to the 
production of heavy mesons in e + e~ or pp collisions. 

The last ingredient are semileptonic decays of heavy mesons. Only re- 
cently the CLEO [3] and BABAR |4j collaborations has measured precise 
spectra of electrons/positrons coming from the decays of D and B mesons. 
This is done by producing specific resonances: ^(3770) which decays into 
D and D mesons (CLEO) and T(4S*) which decays into B and B mesons 
(BABAR). In both cases the heavy mesons are almost at rest, so in practice 
one measures the meson rest frame distributions of electrons/positrons. 

In this presentation the results for production of electrons have been 
obtained within the /^-factorization approach. At relatively low RHIC en- 
ergies intermediate x-values of gluon distributions become relevant. The 
Kwiecinski unintegrated gluon (parton) distributions seem the best suited 
in this context [3J. We use both Peterson [U] and so-called perturbative 
|10j fragmentation functions. The electron/positron decay functions fitted 
recently [12] to the recent CLEO and BABAR data are used. 

Recently the PHENIX collaboration has measured dilepton invariant 
mass spectrum from to 8 GeV in pp collisions at y/s = 200 GeV [7J. Up 
to now, production of open charm and bottom was studied only in inclusive 
measurements of charmed mesons [|5j and electrons [6] and only inclusive 
observables were calculated in pQCD approach [TT| [T2]. 

Some time ago we have studied kinematical correlations of cc quarks [2] , 
which is, however, difficult to study experimentally. High luminosity and in 
a consequence better statistics at present colliders opens a new possibility to 
study not only inclusive distributions but also correlations between outgoing 
particles. Kinematical correlations constitute an alternative method to pin 
down the cross section for charm and bottom production. 

In this presentation I shall show some selected results obtained in |12l 
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Recently we have studied a new mechanism of exclusive production of 
cc pairs [20] . In such a process a single pair of cc is produced together with 
associated two protons. We shall comment on a possibility to identify the 
peculiar mechanism. The original presentation at the conference included 
also inclusive diffractive processes. 

With growing energy the heavy quark-antiquark production becomes 
sensitive to lower- x gluon distributions. At high energies a possibility of 
two pairs of cc production opens up |21j . We comment on which areas of 
the phase space are potentially interesting in order to pin down the double 
parton scattering contribution. 



2.1. Formalism 

We consider the reaction h% + h% — »■ Q + Q + X, where Q and Q are 
heavy quark and heavy antiquark, respectively. 



Fig. 1. Basic diagrams relevant for gluon-gluon fusion (left panel) and quark- 
antiquark annihilation (right panel) in the ^-factorization approach. 

In the fet-factorization approach the differential cross section reads: 



where Ti(x%, k\ t ) and J 7 j(x2, k| t) are the so-called unintegrated gluon (par- 
ton) distributions. Leading-order matrix elements for off-shell gluons |15] 
were used. The two-dimensional Dirac delta function assures momentum 



2. Inclusive production of cc 
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conservation. The unintegrated parton distributions are evaluated at: 

mi t / \ mi± i \ mi t / \ ni2,t / \ 

x i = ~t ex P(yi) + —rf- exp(y 2 ), x 2 = exp(-yi) + exp(-y 2 ), 

where = + m^. 

Introducing new variables: = K\ t t + /t 2 ,t ; 
= H\,t — K2,t one can write: 

dan 



dyidy2d 2 pijd 2 p2,t J ^ * 47r 2 167t 2 (xix 2 s) 2 ^ 8 - 7 ' 

Fi{xi,K 2 l t ) Fj{x2,nl t ) . (2) 

This formula is very useful to study correlations between the produced heavy 
quark Q and heavy antiquark Q [II] . 

At the Tevatron and LHC energies the contribution of the gg —> QQ 
subrocess is more than one order of magnitude larger than its counterpart 
for the qq — > QQ subprocess. At RHIC energy the relative contribution 
of the qq annihilation is somewhat bigger. Therefore in the following we 
shall take into account not only gluon-gluon fusion process but also the 
quark-antiquark annihilation mechanism. 

The production of electrons/positrons is a multi-step process. The whole 
procedure of electron/positron production can be written in the following 
schematic way: 

da e da Q 

Tl2~ = Tl2~ D Q^D ® f D ^e , (3) 
dyd A p dyd z p 

where the symbol ® denotes a convolution of the different distributions. 
The first term is responsible for production of heavy quarks/antiquarks. 
Next step is the process of formation of heavy mesons. We follow a phe- 
nomenological approach and take e.g. Peterson [9] and Braaten et al. [10] 
fragmentation functions with parameters from the literature [16] . The elec- 
tron decay function accounts for the proper branching fractions. 

The inclusive distributions of hadrons can be calculated as a convolution 
of inclusive distributions of heavy quarks/antiquarks and Q — > h fragmen- 
tation functions: 

dcr(yi,Pi t ,y 2 ,P2tA) _ f Dq^ h (z{) Dq^ r {z 2 ) du{y 1 ,p%y2,V%A) dz dz 



dyidp^ t dy 2 dp% t d(t) J z x z 2 dyidp^ t dy 2 dp^ t d(t) 

(4) 

where: pf t = p^t = with meson longitudinal fractions z\,Z2 £ (0, 1). 

We use decay functions fitted recently [12] to the CLEO and BABAR 
data. In our approach the electrons (positrons) are generated isotropically 
in the heavy meson rest frame. 
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2.2. Results 

2.2.1. Single electron spectra 

Before we start presenting our single particle distributions let us focus 
on the decay functions discussed shortly above. In Figj2]we show our fit [12] 
to the CLEO and BABAR data. The good quality fit of the data allows to 
obtain reliable predictions for electron/positron single particle spectra. 




P. (GeV) p. (GeV) 



Fig. 2. Our fit to the CLEO [3] and BABAR g] data. 

Now we shall focus on transverse momentum distribution of electrons/positrons 
measured recently by the STAR and PHENIX collaborations at RHIC [TJ [2] . 
In Fig|3j as an example, we show results obtained with the Kwiecihski 
UPDFs [8j. In Ref. [12] we have discussed in addition other UGDFs. Uncer- 
tainties due to different combinations of factorization and renormalization 
scales as well as due to different choices fragmentation functions are shown 
in FigSl In these calculations we have included both gluon-gluon fusion 
as well as quark-antiquark annihilation. In the last case we use matrix 
elements with on-shell formula but for off-shell kinematics (the discussion 
of this point can be found in our earlier paper [H]). In Ref. [12] we have 
discussed also uncertainties due to the choice of quark masses. 

Study of nonphotonic and hadron correlations allows to "extract" 
a fractional contribution of the bottom mesons B/(D + B) as a function 
of electron/positron transverse momentum [T7]. Recently the STAR col- 
laboration has extended the measurement of the relative B contribution to 
electron/positron transverse momenta ~ 10 GeV [18] . In Figj5] (Kwiecihski 
UPDFs) shown are results for different scales and different fragmentation 
functions. There is a strong dependence on the factorization and renor- 
malization scale. A slightly better agreement is obtained with the Peterson 
fragmentation functions. 



G 
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P, (GeV) p t (GeV) 



Fig. 3. Transverse momentum distribution of electrons/positrons obtained with 
the Kwiecihski UPDFs. We show separately contributions of the gluon-gluon fusion 
(black) and quark-antiquark annihilation (grey). On the left side results with the 
Peterson fragmentation functions and on the right side with BCFY fragmentation 
functions. 

2.2.2. Electron-positron correlations 

When calculating correlation observables we have included also processes 
shown in FigJU] and FigJTl The photon-photon induced processes were first 
included in Ref. [13] . The central exclusive diffractive process shown in 
FigJ7]was proposed in Ref. [20]. 

In FigJH]we show e + e~ invariant mass distributions calculated with the 
Kwiecinski (left) and KMR (right) UGDFs. One can clearly see that both 
the Kwiecinski and KMR [22] UGDFs give fairly good description of the 
data for M e + e - > 3 GeV. At small invariant masses the Kwiecinski UGDF 
underestimates the PHENIX data and the KMR UGDF starts to overesti- 
mate the data points below M e + e - = 2 GeV. 

In Figj9] we show uncertainties related to the contribution of semilep- 
tonic decays. The left panel presents uncertainties due to the factorization 
scale variation as described in the figure caption. The right panel shows 
uncertainties due to the modification of the heavy quark masses. 

If the detector can measure both transverse momenta of electron/positron 
and their directions, one can construct a distribution in transverse momen- 
tum of the dielectron pair: pt )S um = Pit + P2t- Our predictions including the 
semileptonic decays and Drell-Yan processes are shown in the left panel of 
FigJIUl Both processes give rather similar distributions. The distributions 
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Fig. 4. Factorization and renormalization uncertainty band of our A^-factorization 
calculation with unintegrated Kwieciriski gluon, quark and antiquark distributions 
for the Peterson fragmentation function (left panel) and BCFY fragmentation func- 
tion (right panel). The open triangles represent the PHENIX collaboration data 
and the solid circles the STAR collaboration data. 




P. (GeV) p, (GeV) 



Fig. 5. The fraction of the B decays for the Kwiecihski UPDFs. The uncertainty 
band due to the choice of the scales is shown for the Peterson (left) and Braaten 
et al. (right) fragmentation functions. Both gluon-gluon fusion as well as quark- 
antiquark annihilation are included in this calculation. 
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Fig. 6. Diagrammatic representation of processes initiated by photon-photon sub- 
processes: double-elastic, double-inelastic, inelastic-elastic and elastic-inelastic. 




Fig. 7. The mechanism of exclusive double-diffractive production of open charm. 




Fig. 8. Dielectron invariant mass distribution for pp collisions at t/s = 200 GeV for 
the Kwiecihski (left) and KMR (right) UGDFs. Different contributions are shown 
separately: semileptonic decay of charm by the blue solid line, semileptonic decay of 
bottom by the red solid line, Drell-Yan mechanism by the long dashed line, gamma- 
gamma processes by the blue dashed line and the central diffractivc contribution 
by the green dotted line. In this calculation we have included azimuthal angle 
acceptance of the PHENIX detector [7]. 



of this type were not measured so far experimentally. The distribution in 
Pt,sum is not only a consequence of gluon transverse momenta but invlolves 
also fragmentation process and semileptonic decays. With good azimuthal 
resolution of detectors one could also construct distribution in azimuthal 
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M„ (GeV) 




Fig. 9. The uncertainties of theoretical calculations. The left panel shows the 
factorization scale uncertainties, the lower curve corresponds to = m\ t + 

m\ t and the upper curve to n 2 R = kf, \j? F = ^itiq, where k t is gluon transverse 
momentum. The right panel shows the quark mass uncertainties as indicated in 
the hgure. 



angle between electron and positron. Corresponding predictions are shown 
in the right panel of FigllOl One can see an interesting dependence on the 
invariant mass of the dielectron pair - the smaller the invariant mass the 
large the decor relation in azimuthal angle. 




Fig. 10. Distribution in transverse momentum of the dielectron pair (left) and 
in azimuthal angle between electron and positron (right) for semileptonic decays 
(solid line) and Drell-Yan processes (dashed line). Here Kwiecinski UGDF and 
Peterson fragmentation function were used. 



3. Exclusive diffractive production of cc 

Central exclusive mechanisms of cc production at high energies shown 
in FigJT] constitutes a special category of diffractive processes. In this case 
the central system X is produced in the color singlet state. This leads to 
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rapidity gaps between forward/backward produced protons and the central 
system. The QCD mechanism of central exclusive heavy quark-atiquark 
dijets (in particularly bb) is a source of the irreducible background to the 
exclusive Higgs boson production [231 123] - Central exclusive production of 
cc and bb pairs was studied in detail in our previous papers [201 1231 121] • 
In these calculations the pp — > p(qq)p reaction was considered as a 4-body 
process with exact kinematics. The applied perturbative model of theoreti- 
cal predictions is based on the Khoze-Martin-Ryskin (KMR) approach used 
previously for the exclusive Higgs boson production [25J. Total cross sec- 
tions and differential distributions for heavy quarks are calculated by using 
fct-factorization approach with the help of off-diagonal unintegrated gluon 
distribution functions. 

This QCD model works very well in the case of exclusive dijets and 
charmonia production (26J EH EEl ES] ■ Here we discuss the production of 
cc pairs. In practice, however, one measures rather charmed mesons. The 
measurement and its interpretation is therefore more complicated and will 
be not discussed here. Such experimental analyses are being performed 
now at the Tevatron and could be also available in Run II at RHIC. In this 
context it is very interesting to compare the mechanism of central exclu- 
sive production of charm quarks with standard inclusive single and central 
diffractive processes. 

As in the KMR approach [251 [30] the amplitude of the exclusive central 
diffractive qq pair production pp — > p{qq)p can be written as 

^v, = | • 0^ 9 / d\ u v^( qi , q2 , h,k 2 ) 

ffi( x i,x'i,Qlt,qlt, h)f° S 2 (x 2 ,x' 2 , gj t , gj t , t 2 ) 

X = o 2 ' V 3 ) 

lit Qi,t lit 

where X q , X q are helicities of heavy q and q, respectively, t\ }2 are the momen- 
tum transfers along each proton line, qij,q 2 ,t,xi,2 an d qo,t, x'\ ~ x 2 <C x\ t2 
are the transverse momenta and the longitudinal momentum fractions for 
active and screening gluons, respectively. Above fg\/ 2 are t ne off-diagonal 
UGDFs related to both nucleons. The vertex factor V£ 1 ^(qi,q 2 ,k%,k 2 ) is 
the production amplitude of a pair of massive quark q and antiquark q with 
helicities \ q , Xq and momenta k±, k 2 , respectively. The longitudinal momen- 
tum fractions of active gluons are calculated based on kinematical variables 
of outgoing quark and antiquark: x% = exp(+?/q) + exp(+yq) and 

x 2 = ^^exp(— y q ) + ^^exp(— y q ), where m q j and m ? - f are transverse 
masses of the quark and antiquark, respectively, and y q and y q are corre- 
sponding rapidities. 
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The off-diagonal UGDFs can be approximated as [31] 

ff(x', xi )2 , 91,24, qlp H%) ^ Rg fg{x\,2, q1,2U /4") • (6) 

The factor i? 9 here cannot be calculated from first principles in the most 
general case of off-diagonal UGDFs. It can be estimated only in the case of 
off-diagonal collinear PDFs when x' <C x and xg = x _A (l — x) n and then 
Rg = r r\\+ff ■ ^ n ^ ne considered kinematics the diagonal unintegrated 
densities can be written in terms of the conventional (integrated) densities 
xg(x,qf) as [M] 



/ 9 (*,^V) = Q^[xg(x,ig)yjT g (<g,ij?)] , (7) 

where T g is the conventional Sudakov survival factor which suppresses real 
emissions from the active gluon during the evolution. 

The hard subprocess g*g* — > qq amplitude V^ 2 _(qi,q2,ki,k2) reads 



(qi — ki) 2 — m 2 (q\ — A^) 2 — m 2 



(8) 



In the present calculations we use the GJR08 set of collinear gluon dis- 
tributions [33]. In the analogy to the CEP of Higgs boson, where renor- 
malization and factorization scales /i 2 = /i 2 ^ = fj, F = Mfj are preferred, we 
take [i 2 = M 2 S . Absorption corrections to the bare pp — > p(qq)p amplitude, 
which are necessary to be taken into account (to ensure exclusivity of the 
process) , are included approximately by multiplying the cross section by the 
gap survival factors Sq = 0.1 for RHIC and Sq = 0.03 for the LHC energy. 
More details about exclusive production of heavy quarks can be found in 
our original paper |20j . Let us come now to presentation of our results. 

In Fig. [11] we show rapidity distribution of c quarks from the exclusive 
mechanism (solid lines) shown already in Fig. 0. We show the results for LO 
(upper curves) and NLO (lower curves) collinear gluon distributions [33j . 
We observe large difference of results for LO and NLO gluon distributions 
especially at LHC. For comparison we show the contribution of inclusive 
central diffractive component discussed in detail in [19]. In this calculation 
we have included gap survival factors Sq = 0.1 for y/s = 500 GeV and Sg 
= 0.03 for \/s = 14 TeV. The cross section for the exclusive mechanism is 
similar to that for the inclusive central diffractive mechanism. The exclusive 
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production starts to dominate only at large c quark rapidities. Therefore a 
measurement of the cross section with double (both side) rapidity gaps may 
be not sufficient to single out the exclusive mechanism. Clearly other cuts 
would be necessary. 




Fig. 11. Distributions in rapidity of c quark/ antiquark for the exclusive component 
at ■</£ = 500 GeV (left panel) and <fs = 14 TeV (right panel) . G JR08 collinear gluon 
distributions were used to obtain the unintegrated gluon distribution according to 
the KMR prescription. For comparison we show the inclusive central diffractive 
contribution (dashed line). 



Distributions in the c quark (c antiquark) transverse momentum are 
shown in Fig. [12j At RHIC energy distributions for both mechanisms have 
very similar shape. However, at LHC nominal energy we observe that inclu- 
sive central diffractive component extends to higher transverse momentum 
than that for the exclusive central diffractive one. In order to identify the ex- 
clusive component a much more precise analysis of kinematical correlations 
between quark and antiquark is needed. A detailed Monte Carlo studies of 
final states of both mechanisms could help to find a criterion to separate 
experimentally the two dynamically different components. 
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Pt (GeV) p t (GeV) 



Fig. 12. Distributions in transverse momentum of c quark/ antiquark for the ex- 
clusive component at s/s = 500 GeV (left panel) and s/s = 14 TeV (right panel). 
GJR08 collinear gluon distributions were used to obtain the unintegrated gluon 
distribution. For comparison we show the inclusive central diffractive contribution 
(dashed line). 



4. Double parton scattering production of cccc 

4-1. Framework 

The double-parton scattering has been recognized long ago. Several 
estimates of the cross section for different processes have been presented in 
recent years. In the present analysis we discuss production of (cc)(cc) four- 
parton final state which has not been discussed so far but is particularly 
interesting especially in the context of experiments being carried out at 
LHC. 

The double-parton scattering formalism proposed so far assumes two 
single-parton scatterings. In a simple probabilistic picture the cross section 
for double-parton scattering can be written as: 

a DPS (pp -)• ccccX) = —^—<T SPS (pp -> ccX x ) ■ a SPS (pp -»• ccX 2 ). (9) 

2(J eff 

This formula assumes that the two subprocesses are not correlated and 
do not interfere. At low energies one has to include parton momentum 
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conservation i.e. extra limitations: x\ + x 3 < 1 and x 2 + £4 < 1, where 
x\ and X3 are longitudinal momentum fractions of gluons emitted from one 
proton and X2 and X4 their counterpairs for gluons emitted from the second 
proton. Experimental data provide an estimate of c e // i n the denominator 
of formula @. In our analysis we take <r e // = 15 mb. 




Fig. 13. Mechanism of (cc)(cc) production via double-parton scattering. 

The simple formula ([9|) can be generalized to include differential distri- 
butions. In the same approximation differential distribution can be written 
as 

da 1 da da 

dyidy 2 d 2 pitdy 3 dy 4l d 2 p2t 2a eff dyidy 2 d 2 p lt dy 3 dyA t d 2 p 2t ' 

which reproduces formula ([9|). This cross section is formally differential in 
8 dimensions but can be easily reduced to 7 dimensions noting that physics 
of unpolarized scattering cannot depend on azimuthal angle of the pair or 
on azimuthal angle of one of the produced c (c) quark (antiquark). The 
differential distributions for each single scattering step can be written in 
terms of collinear gluon distributions with longitudinal momentum fractions 
xi, X2, xs and x^ expressed in terms of rapidities yi, y 2 , 2/3, 2/4 and transverse 
momenta of quark (or antiquark) for each step. 

A slightly more general formula for the cross section can be written 
formally in terms of double-parton distributions (dPDF), e.g. F gg , F qq , etc. 
In the case of heavy quark (antiquark) production at high energies: 

da DPS = — — FggiX!^^, f4)F gg (x 2 ,X 4 , f4) 

1(7 • If 

da gg ^ C c(x 1 ,x 2 , nl)dagg^ C c(x 3 , 24, fj%) dx 1 dx 2 dx 3 dx 4 . (11) 

The double-parton distributions in Eq. (jlip are not well known. Usually 
one assumes the factorized form and expresses them via standard distribu- 
tions for SPS. Even if factorization is valid at some scale, QCD evolution 
leads to a factorization breaking |21j . 
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In this presentation we shall apply the commonly used factorized model. 
Some refinements are presented in |21| . 



4-2. Results 

In Fig. [T3]we compare cross sections for single cc and DPS cccc produc- 
tion as a function of pp center-of-mass energy. At low energies the cross 
section for cc is much larger. For reference we show the proton-proton total 
cross section as a function of energy. At low energy the cc or cccc cross 
sections are much smaller than the total cross section. At higher ener- 
gies the contributions approach the parametrized total cross section. This 
shows that inclusion of unitarity effect and/or saturation of parton distri- 
butions may be necessary. At LHC energies the cross section for both terms 
becomes comparable. This is a new situation when the double-parton scat- 
tering gives a large contribution to inclusive charm production. This issue 
was not discussed so far in the literature. 




Fig. 14. Total LO cross section for single-parton and double-parton scattering 
as a function of center-of-mass energy (left panel) and uncertainties due to the 
choice of (factorization, renormalization) scales (right panel) . We show in addition 
a parametrization of the total cross section in the left panel. 

As an example in Fig. [15] we present single c (c) distributions. Within 
approximations discussed here the single-parton distributions are identical 
in shape for cc and cccc. This means that double-scattering contribution 
produces naturally an extra energy-dependent if-factor to be contrasted 
with approximately energy-independent /C-factor due to higher-order cor- 
rections. A strong dependence on the factorization and renormalization 
scales can be observed. 

So far we have discussed only single particle spectra of c or c. A bet- 
ter test of DPS could be done by studying correlation observables. The 
correlations between c and c has been studied e.g. in [J3]. In Fig. [Tol we 
show distribution in the difference of c and c rapidities (left panel) as well 
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Fig. 15. Distribution in rapidity (left upper panel) and transverse momentum 
(right upper panel) for different UGDFs and associated ucertainties related to the 
choice of renormalization and factorization scales (lower panels) for core quarks 
at yfs = 7 TeV. 



as in the cc invariant mass M cg (right panel). We show both cases: when cc 
are emitted in the same parton scattering (c\C2 or C3C4) and when they are 
emitted from different parton scatterings [c\C4 or C2C3). We observe a long 
tail for large rapidity difference as well as at large invariant masses of cc. 
Such distributions for quarks and antiquarks cannot be directly measured. 
Instead their counterparts for mesons can be studied. This was discussed 
in more detail in our recent original paper |21j . 

As the last example in Fig. [T7] we present distribution in the transverse 
momentum of the cc pair |pj_ c g|, where p± C 5 = P±l + P±s- F° r compari- 
son this is a Dirac delta function in the leading-order approximation to cc 
production. In contrast, double-parton scattering mechanism gives a broad 
distribution extending to large transverse momenta. NLO corrections obvi- 
ously destroy the <5-like leading-order correlation. Similar distributions for 
DD seem useful observables to identify the DPS contributions [21j . 

So far we have calculated cross section in a simple leading-order ap- 
proach. A better approximation would be to include multiple gluon emis- 
sions. This can be done e.g. in soft gluon resummation or in the kt- 
factorization approach. For example the second approach does not lead 
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• DPS pp->ccccX 

FLO: gg-.cc 
: CTEQS LO 



\ls = 7 TeV 




: DPS pp-»OOOCX 

FLO: gg^cc 




Fig. 16. Distribution in rapidity difference (left panel) and in invariant mass of 
the cc pair (right panel) at y/s — 7 TeV. 



to large changes in shape of neither distributions in rapidity nor of dis- 
tributions in transverse momentum of c (c) (see e.g. [13]) compared to 
the collinear approach. It is expected, however, to change distributions in 
transverse momentum of cc or in azimuthal angle between c and c |14j . 



i DPS pp- 


KCCCX Vs 


= 7 TeV : 


!T LO: gg^cl? 






it I/JS8.0 


c,n 4 


or c 2 C 3 n 






or c 3 c a 






CTE06LO i 









2 4 6 8 10 12 14 16 18 20 



Fig. 17. Distribution in transverse momentum of cc pairs from the same parton 
scattering and from different parton scatterings at y/s — 7 TeV. 



5. Conclusions 

We have calculated inclusive spectra of nonphotonic electrons/positrons 
for RHIC energy in the framework of the fcj-factorization. We have con- 
centrated on the dominant gluon-gluon fusion mechanism and used unin- 
tegrated gluon distribution functions from the literature. Special emphasis 
was devoted to the Kwiecinski unintegrated gluon (parton) distributions. 
In this formalism, using unintegrated quark and antiquark distributions, 
one can calculate also the quark-antiquark annihilation process including 
transverse momenta of initial quarks/antiquarks. 
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When calculating spectra of charmed (D, D*) and bottom (B, B*) 
mesons we have used Peterson and Braaten et al. fragmentation functions. 
We have used recent fits to the CLEO and BABAR collaborations data for 
decay functions of heavy mesons. 

Our results have been compared with experimental data measured re- 
cently by the PHENIX and STAR collaborations at RHIC. A reasonable 
description of the data at large transverse momenta of electrons/positrons 
has been achieved. We have discussed uncertainties related to the choice of 
the factorization and renormalization scales as well as those related to the 
fragmentation process. Although the uncertainty bands are rather large, 
there seems to be a missing strength at lower electron/positron transverse 
momenta. 

We have discussed also correlations of charmed mesons and dielectrons 
at the energy of recent RHIC experiments. We have calculated the spec- 
tra in dielectron invariant mass, in azimuthal angle between electron and 
positron as well as the distribution in transverse momentum of the pair. 
The uncertainties due to the choice of UGDFs, choice of the factorization 
and renormalization scales, have been discussed. We have obtained good 
description of the dielectron invariant mass distribution measured recently 
by the PHENIX collaboration at RHIC. 

At RHIC the contribution of electrons from Drell-Yan processes is only 
slightly smaller than that from the semileptonic decays. The distributions 
in azimuthal angle between electron and positron and in the transverse 
momentum of the dielectron pair from both processes are rather similar. 
It was found that the distribution in azimuthal angle strongly depends on 
dielectron invariant mass. 

We have also included exclusive central-diffractive contribution discussed 
recently in the literature. At the rather low RHIC energy it gives, however, 
a very small contribution to the cross section and can be safely ignored. 

The QED double-elastic, double-inelastic, elastic-inelastic and inelastic- 
elastic processes give individually rather small contribution but when added 
together are not negligible especially at low dielectron invariant masses 
where some strength is missing. 

The exclusive production of cc pairs is interesting by itself. We have 
discussed corresponding formalism as well as some results for RHIC and 
LHC energies. However, experimental identification of the mechanism may 
be not easy as the final hadronic state is more complicated and will compete 
with inclusive central diffractive production of cc. 

We have discussed also production of two pairs of cc. We have found 
very quickly rising cross section for the two-pair production as a function 
of center-of-mass energy. The two-pair production must therefore give a 
sizeable contribution to inclusive charm production. This point requires 
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further studies. 

We have discussed some promissing observables which seem useful in 
identifying the DPS production of two pairs of cc. In Ref.[21| we have 
considered also corresponding observables for charmed mesons. Another 
option would be to study production of the same-sign charged leptons. We 
expect that semileptonic decays are the main source for semi-hard muons 
or electrons. Furthermore this contribution can, in principle, be separated 
experimentally by taking into account that the secondary vertices are shifted 
with respect to the primary ones. This should allow a separation of the 
semileptonic "signal" from other possible sources of dilepton continuum. 
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